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Abstract 

We reconsider the effects of nnstable particles on the prodnction and destruction of the 
primordial light elements, with a view to reconciling the high primordial ^Li abundance 
deduced from Big Bang Nucleosynthesis (BBN), as implied by the baryon-to-photon ratio 
now inferred from the anisotropies of the Cosmic Microwave Background (CMB), with the 
lower abundance of ^Li observed in halo stars. The potential destruction of ^Li is strongly 
constrained by observations of Deuterium (D), ^He and ®Li. We identify ranges for the 
unstable particle abundance and lifetime which would deplete ^Li while remaining consistent 
with the abundance of ®Li. However, in these regions either the D abundance is unacceptably 
low or the ratio ^He/D is unacceptably large. We conclude that late particle decay is unable 
to explain both the discrepancy of the calculated ^Li abundance and the observed ^Li plateau. 
In the context of supersymmetric theories with neutralino or gravitino dark matter, we 
display the corresponding light-element constraints on the model parameters. 
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1 Introduction 


The observed abundances of light elements are generally in good agreement with the predic¬ 
tions of Big Bang Nucleosynthesis (BBN) calculated assuming a homogeneous Robertson- 
Walker-Friedman cosmology [1]. Within this framework, the snccess of BBN calculations 
imposes important constraints on the nnmber of light particle species and on the baryon-to- 
photon ratio f] [2], Recent high-precision measnrements of the cosmic microwave backgronnd 
(CMB) radiation by WMAP and other experiments now complement the BBN in important 
ways [3]. For example, they impose strong constraints on r] and weaker constraints on the 
nnmber of light particle species. Of particular interest are the very precise predictions of the 
light element abnndances from BBN that are made with the CMB valne of r] [4]. 

Based on these predictions, there is now tension between some observed light-element 
abnndances and those that wonld be calculated using the CMB valne of r] and assnming no 
additional light particles beyond three light nentrino species. In particnlar, the prediction for 
the primordial abnndance of ’^Li made using the CMB valne of r] [5, 6, 7, 8, 9] is somewhat 
higher (by a factor of 2 - 3) than the primordial abnndance inferred from astrophysical 
observations. The significance of this discrepancy should not be over-emphasized, in view of 
the potential systematic errors in the interpretation of the astrophysical data. However, it has 
stimnlated theoretical explorations of mechanisms for modifying the CMB/BBN prediction 
for ^Li, for example via the late decays of massive particles [10, 11]. 

The effects of snch decays have been studied extensively, and the constraints imposed on 
them by the observed abnndances of ^He, Denterium and ®Li are well nnderstood [12, 13, 14]. 
Previous studies had shown no incompatibilities between these constraints and the snggestion 
that late-decaying particles might have modified the BBN prediction for the abnndance of 
^Li so as to agree better with the astrophysical observations. Snpersymmetric models with 
conserved R parity natnrally predict such a particle, either the gravitino if it is not the 
lightest snpersymmetric particle (LSP), or a nentralino or stan slepton if the gravitino is the 
LSP [15, 16], which may well have the appropriate abnndance and lifetime to affect both 
the ®Li and ^Li abundances. Later in this paper, we identify these regions in the parameter 
space of a constrained version of the minimal snpersymetric extension of the Standard Model 
(CMSSM), as well as in very constrained versions of the model [17] motivated by snpergravity 
considerations, which predict a gravitino LSP in parts of the parameter space. 

However, we also point ont in this paper that snch scenarios typically yield an abnndance 
of ^He that may be more than an order of magnitnde larger than the Denterium abundance. 
Since it appears that the abnndance of ^He has remained relatively constant in time, whereas 
Denterinm (D) would have been destroyed in stars, it seems unlikely that the ratio of the ^He 
to D abnndances could have been significantly larger than it present valne in the early history 
of the Universe [14]. Imposing this constraint, we find a farther restriction on the parameter 
space of late-decaying heavy-particle models, which exclndes the region of parameter space 
where they could have the desired impact on the primordial abnndance of ^Li. We display 
this and other constraints on the lifetime and abnndance of any massive unstable relic particle 
such as a gravitino, pointing out also the potential impact of weakening the lower limit on the 
primordial Denterinm abnndance. We also display the effects of the ^He constraint on the 
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parameter spaces of models in which heavier supersymmetric sparticles decay into gravitinos. 
The overall reductions in the allowed parameter spaces are often small, but they do exclude 
the regions where the abundance of ^Li could be brought into line with the astrophysical 
observations. This may motivate a reassessment of their interpretation. 

2 Is there a problem with Lithium-7? 

The most direct and accurate estimate of the baryon-to-photon ratio rj is currently provided 
by the acoustic structures in the CMB perturbations [3], namely 

7] = 6.14 ±0.25 X 10"^°. (1) 

This range may be used an input into homogeneous BBN calculations [5, 6, 7, 8, 9], yielding 
the following abundances for the elements of principal interest, which are taken from [9]: 

0.2485 ±0.0005, 

2.55t°;2i X 10-^ 

1.01 ± 0.07 X 10■^ 

1.3±0.1 X 10-^^ (2) 

where Yp is the ^He mass fraction, and the other abundances are expressed in terms of their 
numbers relative to H, as shown. 

These abundances may then be compared with the abundances of the same elements that 
are inferred from those observed in the most primitive astrophysical sites [18, 19, 20, 21, 22]: 

0.232 to 0.258, 

2.78 ±0.29 X 10■^ 

1.5 ±0.5 X 10■^ 

1.23l°;i X 10-'°, 

Qtl X 10-^2. (3) 

Note that, for ^He, we have at our disposal only local data coming from star-forming HIT 
regions in the galactic disk [23] or from the proto-solar value [24]. Comparing the two sets 
of abundances, (2) and (3), we see no signihcant discrepancies, except in the case of ^Li and 
®Li. However, whereas ^Li has to be destroyed, one must produce a factor 1000 more ®Li. 
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The value quoted above for the ^Li abundance assumes that Li depletion is negligible in 
the stars observed. Indeed, standard stellar evolution models predict Li depletion factors 
which are very small (less than 0.05 dex) in very metal-poor turnoff stars [25]. However, 
there is no reason to believe that such simple models incorporate all effects which could lead 
to depletion, such as rotationally-induced mixing and/or diffusion. Including these effects, 
current estimates for possible depletion factors are in the range ~ 0.2-0.4 dex [26]. However, 
the data sample [27] used in deriving the abundance in (3) shows a negligible intrinsic spread 
in Li, leading to the conclusion that depletion in these stars is in fact quite low (less than 
0.1 dex). 

Another important source of potential systematic uncertainty is related to the assumed 
surface temperature of the star. A recent study [28] found significantly and systematically 
higher temperatures for stars used in ^Li observations, specihcally at low metallicity. This 
result leads to ^Li/H = (2.34 ± 0.32) x 10“^'^, which is still, however, nearly a factor of 2 
smaller than the BBN/WMAP prediction. We note, hnally, that another potential source 
for theoretical uncertainty lies in the BBN calculation of the ^Li abundance. However, this 
too has been shown to be incapable of resolving the ^Li discrepancy [6, 29]. 

The comparison of (2) and (3) also shows a discrepancy for ®Li. It has generally been 
assnmed that ®Li is prodnced in post-BBN processes snch as galactic cosmic-ray nucleosyn¬ 
thesis [30]. Until recently, the abundance of ®Li had been observed only in a few metal-poor 
halo stars with metallicity [Fe/H] larger than -2.3. However, new observations of this isotope 
have now been obtained in halo stars. New valnes of the ratio ®Li/’^Li have been measured 
with UVES at the VLT-UT2 Kueyen ESO telescope, in halo stars with metallicity ranging 
from -2.7 to -0.5 [21, 22]. These observations indicate the presence of a plateau in ®Li/H 
~ 10“^^, suggesting a pregalactic origin for the formation of ®Li [31]. These data provide 
interesting new constraints in the present context. 

The relation of the observed ^Li abundance to its primordial value may be debated, but for 
the moment we take the discrepancy between the observed and calculated ^Li abundances 
at face valne, and explore its possible theoretical interpretation in terms of late-decaying 
massive particles, depleting ^Li and possibly producing ®Li, withont negative effects on D 
and ^He. 

3 Possible Impact of Late-Decaying Particles 

The decays of massive particles X with lifetimes > 10^ s could, in principle, have modihed 
the BBN predictions in either of two ways. First, their decay products would have increased 
the entropy in the primordial plasma, implying that its value during BBN was lower than 
that inferred from the CMB. However, this effect is negligible compared with the second 
effect, which is the modihcation of the BBN light-element abnndances by the interactions of 
decay products [12]. 

The latter possibility has been explored recently as a mechanism for reducing the pri¬ 
mordial ^Li abnndance [10, 11]. The existence of such late-decaying massive particles X 
is a generic possibility in supersymmetric models, in particular [15, 16]. Examples include 
a massive gravitino weighing ~ 100 GeV, if it is not the lightest supersymmetric particle 
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(LSP), or some other next-to-lightest supersymmetric particle (NSP) if the gravitino is the 
LSP. Cosmological constraints on such scenarios have been explored previously [12, 13, 14], 
However, the potential ^Li problem motivates a re-examination of the astrophysical and cos¬ 
mological constraints on such scenarios. In particular, we wish to determine whether or not 
a possible solution to this problem can be found in the context of motivated and well-studied 
supersymmetric models. 

3.1 The D and ^He Abundances 

Fig. 1 shows the current constraints in the (tx, Cx) plane, updating those shown in [13]. The 
green lines are the contours 

(1.3or 2.2) X 10-® < ^ < 5.3 x 10“^ (4) 

H 

The first of the lower bounds is the higher line to the left of the cleft, and represents the 
very conservative lower limit on D/H assumed in [13]. The range 1.3 - 5.3 xl0“^ effectively 
brackets all recent observations of D/H in quasar absorption systems. The second of the 
lower bounds is the lower line on the left side, and represents what now seems a reasonable 
lower bound, which is obtained from the 2-a lower limit in (3). The upper bound in (4) is 
the line to the right of the cleft, and is the same upper limit as was used in [13]. A priori, 
there is also a narrow strip at larger C,x and tx where the D/H ratio also falls within the 
range (4), but this is excluded by the observed "^He abundance. 

The solid red lines in the npper right part of Fig. 1 correspond to the limits 

Yp > 0.227 or 0.232, (5) 

where the lower nnmber (corresponding to the higher line) was used in [13], and the higher 
nnmber (corresponding to the lower line) is a lower limit that has been advocated recently 
[18]. It is apparent that, for onr pnrposes, the third signihcant hgure in the "‘He abnndance 
is nnimportant: the narrow D/H strip is in any case exclnded, and there are always stronger 
bonnds on C,x at large tx- 

3.2 The ®Li Abundance 

As said above, recent observations of ®Li in halo stars have provided new insight into the 
origin and the evolntion of this isotope. We recall that ®Li is a pnre prodnct of spallation, 
and many stndies have followed the evolntion of ®Li in onr Galaxy [30]. Of particular 
importance in this context is the a + a reaction that leads to the synthesis of this isotope 
as well as ^Li, and is efficient very early in the evolntionary history of the Galaxy. The 
new valnes of ®Li/^Li that have been measnred in halo stars with UVES at the VLT-UT2 
Kneyen ESO telescope indicate the presence of a platean in ®Li, which snggests a pregalactic 
origin for the formation of this isotope. The evolntion of ®Li with redshift was calculated 
[31] following an initial bnrst of cosmological cosmic rays np to the formation of the Galaxy. 
This process is able to prodnce the reqnired abnndance of ®Li observed in metal-poor halo 
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Figure 1: The constraints imposed by the astophysical observations of^He (red lines), D/H 
(green lines), ^Li (yellow line), ^Li/^Li (blue lines), Li (blue band) and ^He (black lines). 
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stars without the additional over-production of ^Li. In this paper, we have to consider the 
new constraint brought by the existence of this plateau. The potential destruction of ^Li by 
unstable particles must not lead to over-production of ®Li. 

The constraint imposed by the ®Li abundance is shown as a solid yellow line in Fig. 1, 
which is the same as that discussed in [13]. Also shown, as solid blue lines, are two contours 
representing possible npper limits on the ®Li/^Li ratio: 

®Li 

— < 0.07 or 0.15, (6) 

‘ ijl 

with the npper (lower) contonr corresponding to the npper (lower) nnmber in (6). The lower 
nnmber was nsed in [13] and represented the npper limit available at the time, which was 
essentially based on mnltiple observations of a single star. The most recent data [21, 22] 
inclndes observations of several stars. The Li isotope ratio for most metal-poor stars in the 
sample is as high as 0.15, and we display that npper limit here. The main new effect of this 
constraint is to disallow a region in the near-vertical cleft between the npper and lower limits 
on D/H, as seen in Fig. 1. 

3.3 The ^Li Abundance 

The main region of interest in Fig. 1 is the bine shaded band that represents the inferred ^Li 
abnndance: 

7Li 

0.9 X 10“^° < — < (2 or 3) X 10"^°, (7) 

H 

with the ^Li abnndance decreasing as (x increases and the intensity of the shading changing 
at the intermediate valne. In [13], only the lower bonnd was used due the existing discrep¬ 
ancy between the primordial and observationally determined valnes. It is apparent that ^Li 
abnndances in the lower part of the range (7) are possible only high in the Denterinm cleft, 
and even then only if one uses the recent and more relaxed limit on the ®Li/^Li ratio (6). 
Values of the ^Li abundance in the npper part of the range (7) are possible, however, even 
if one uses the more stringent constraint on ®Li/^Li. In this case, the allowed region of 
parameter space would also extend to lower tx, if one conld tolerate valnes of D/H between 

1.3 and 2.2 x 10"^ 

For the convenience of the snbseqnent discussion, the region of the (tx,(x) plane that 
is of interest for lowering the ^Li abnndance is shown alone in panel (a) of Fig. 2. The 
bine region in the arc at low tx is the region exclnded by the stronger lower limit on the 
Denterinm abnndance: D/H > 2.2 x 10“^, and the red region to its right is the extra domain 
that is exclnded by the ^He/D ratio, as we discuss below. 

4 The Importance of the ^He Abundance 

We now come to the constraint from ^He, namely that the ^He/D ratio is absnrdly high in 
the Denterinm cleft. Panel (b) of Fig. 2 shows a histogram of the valnes of the ^He/D ratio 
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Figure 2: (a) The region of the {jxXx) plane in which a decaying relic particle could have 
the desired impact on the Li abundance. To derive the blue (darker grey) region, the only 
abundance cuts applied are: 0.9 < ^Li/H x 10^° < 3.0, 1.3 < D/H x 10® < 5.3, and ®Li/^Li 
< 0.15. To obtain the red (lighter grey) region, the lower bound on D/H was increased to 
2.2 X 10“®. (b) A histogram of the ^He/D ratios found in scenarios sampling the region 
displayed in panel (a), with similar colour coding. 


found in a dense sample of scenarios in the interesting regions shown in panel (a). Since 
Deuterium is more fragile than ^He, whose abundance is thought to have remained roughly 
constant since primordial nucleosynthesis when comparing the BBN value to it proto-solar 
abundance, one would expect, in principle, the ^He/D ratio to have been increased by stellar 
processing. Indeed, there is considerable uncertainty in the evolution of ^He [32]. This 
uncertainty is largely associated with the degree to which ^He is produced or destroyed in 
stars. Since D is totally destroyed in stars, the ratio of ^He/D can only increase in time 
or remain constant if ^He is also completely destroyed in stars. The present or proto-solar 
value of ^He/D can therefore be used to set an upper limit on the primordial value. Fig. 1 
displays the upper limits 


^He 

~1T 


< 1 or 2 


( 8 ) 


as solid black lines. Above these contours, the value of ^He/D increases very rapidly, and 
points high in the Deuterium cleft of Fig. 1 have absurdly high values of ^He/D, exceeding 
the limit (8) by an order of magnitude or more. These are the red points producing the 
high-end peak of the histogram shown in panel (b) of Fig. 2, whereas the blue points are 
those excluded by the lower limit D/H > 2.2 x 10“® that is now preferred. We see that these 
points mostly have acceptably low values of ^He/D, though some large values are found near 
the boundary with the red region in panel (a) of Fig. 2. 

It is interesting to note that in the red region which has acceptable D/H and a ^Li/H 
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abundance low enough to match the observed values, the ®Li abundance is relatively high: 
7.3 X 10“^^ <®Li /H< 1.6 X 10“^^. This matches the new ®Li observations quite well, and 
would circumvent the need for an early period of ®Li production by cosmological cosmic rays. 
Unfortunately, ^He/D ranges from 17 - 37 for these parameter values. 

The previous upper limit on rjx [13] corresponded to the constraint mxnx/n^ < 5.0 x 
10“^^ GeV for tx = 10® s. The weaker (stronger) version of the ®He constraint adopted here 
corresponds to 

mx— < 2.0(0.8) X 10-^2 GeV (9) 

for Tx = 10® s. However, the impact of the ®He constraint is even stronger for tx = 10^ s, the 
location of the previous Deuterium cleft. The analysis of [13] would have given mxnx/n^ < 
360 X 10“^^ GeV, whereas the weaker (stronger) ®He constraint adopted here corresponds to 

mx— < 9.3(3.8) x lO’^^ GeV (10) 

for Tx = 10’^ s. 

5 Applications to Supersymmetric Scenarios 

We now discuss some examples of the consequences of the ®He constraint for various super- 
symmetric scenarios in which R parity is conserved. In such models, if the gravitino is not 
the LSP it will generically decay gravitationally with a long lifetime. If the gravitino is the 
LSP, the next-to-lightest supersymmetric particle (NSP) will decay gravitationally into the 
gravitino LSP, again with a long lifetime. 

5.1 Models with an Unstable Gravitino 

We hrst consider the possibility that the gravitino is not the lightest supersymmetric particle 
(LSP), which is instead the lightest neutralino y. In this case, the gravitino is unstable, with 
a lifetime that could well fall within the range considered here. In such a scenario, the light- 
element abundances impose an important upper limit on the possible temperature of the 
Universe, e.g., during reheating after inflation, which we denote by Tr [12, 13, 14]. We recall 
that thermal reactions are estimated to produce an abundance of gravitinos given by [13]: 

!^ = (0.7-2,7)x10-x(^^), (11) 

Assuming that m 3/2 = 100 GeV and tx = 10® s, and imposing the constraints (9), we now 
hnd 

Tr < (0.8 - 2.8) X 10^ GeV, ((0.3 - 1.1) x 10^ GeV) (12) 

for the weaker (stronger) version of the ®He constraint. This becomes an even more signihcant 
constraint on inflationary models, which were already somewhat embarrassed by the previous 
upper limit Tr ~ 2 x 10^ GeV. 
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5.2 Models with Gravitino Dark Matter 

We now consider the possibility that the gravitino is itself the LSP, in which case the next-to- 
lightest supersymmetric particle (NSP) would instead be unstable, also with a long lifetime. 
We will work in the context of the constrained MSSM (CMSSM) [33, 34, 35], so that all 
the scenarios we consider have universal soft supersymmetry-breaking scalar masses mo and 
gaugino masses mi /2 at the GUT scale before renormalization. The magnitude of the Hig- 
gsino mixing parameter |/r| and the pseudoscalar Higgs mass m^ are hxed by the electroweak 
vacuum conditions. We consider scenarios with /i > 0, which are favoured hj Qfj, — 2 and, to 
a lesser extent, b s'j. The scenarios that we study differ in their assumptions about the 
relationship of the gravitino mass m 3/2 to mo, but they all share the common feature that 
the LSP is the gravitino in generic domains of parameter space. 

In the hrst set of scenarios, shown in Fig. 3, we £x the ratio of supersymmetric Higgs 
vacuum expectation values tan (3 = 10, which is among the lower values consistent with our 
hypotheses, and assume (a) m 3/2 = 10 GeV, (b) m 3/2 = 100 GeV, (c) mz /2 = 0.2mo and (d) 
^^ 3/2 = ttiq. In each panel of Fig. 3, we display accelerator, astrophysical and cosmological 
constraints in the corresponding (mi/ 2 ,mo) planes, concentrating on the regions to the right 
of the near-vertical black lines, where the gravitino is the LSP. The vertical black dashed 
and (red) dot-dashed lines represent the lower limits on mi /2 implied for each value of mo by 
the non-observation of a chargino and a Higgs boson at LEP, the latter having a theoretical 
uncertainty 5mi/2 ~ 50 GeV. The (pale green) narrow diagonal strips represent the regions 
where the relic density of the NSP would have lain in the range 0.094 < VLh? < 0.129 favoured 
by WMAP and other measurements of the cold dark matter density, if the gravitino had not 
been the LSP. In fact, the gravitino is always the LSP in the scenarios considered. The NSP 
may be either the lightest neutralino y or the lighter supersymmetric partner of the r lepton: 
this is lighter than the neutralino y below the (red) dotted line. 

Below and to the right of the upper (purple) dashed lines, the density of relic gravitinos 
produced in the decays of other supersymmetric particles is always below the WMAP upper 
limit: VL^i 2 h? < 0.129. To the right of the lower black solid lines, the lifetime of the next- 
to-lightest supersymmetric particle (NSP) falls below 10^ s, and the analysis of [13] cannot 
evaluate the astrophysical constraints from the light-element abundances, in the absence of 
a suitably modified BBN code. The code used in [13], when combined with the observational 
constraints used in [13], yielded the astrophysical constraint represented by the dashed grey- 
green lines in the different panels of Fig. 3. 

These constraints on the GMSSM parameter plane were computed in [16]. For each point 
in the (mi/ 2 ,mo), the relic density of either y or f is computed and (fx is determined using 
Ltxh^ = 3.9 X 10^ GeV Cx- When X = f, (x is reduced by a factor of 0.3, as only 30% of 
stau decays result in electromagnetic showers which affect the element abundances at these 
lifetimes. In addition, at each point, the lifetime of the NSP is computed. Then for each tx, 
the limit on (x is found from the results shown in Fig. 1. The region to the right of this 
curve where r = Cx/Cx™** < 1 is allowed. 

The astrophysical constraints obtained with the newer abundance limits used here yields 
the solid red lines in Fig. 3. The examples where tx and Cx for the NSP decays fall within 
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Figure 3: The (mi/ 2 ,mo) planes for p > 0, tan/3 = 10 and (a) 1713/2 = 10 GeV, (b) 
^^ 3/2 = 100 GeV, (c) m 3/2 = 0.2mo and (d) 1013/2 = itiq. We restrict our attention to the 
regions between the solid black lines, where the gravitino is the LSP and the NSP lifetime 
exceeds 10^ s. In each panel, the near-vertical dashed black (dash-dotted red) line is the 
constraint m^± > 104 GeV (ruh > 114 GeV), the upper (purple) dashed line is the constraint 
kl 3 / 2 h^ < 0.129, and the light green shaded region is that where the NSP would have had 
0.094 < Qh"^ < 0.129 if it had not decayed. The solid red (dashed grey-green) line is the 
region now (previously) allowed by the light-element abundances: r < 1 as described in the 
text. The red (blue) shaded region is that where the abundance could have been improved 
by NSP decays, but which is now excluded by the ^He (D) constraint. 
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the ranges shown in Fig. 2(a), and hence are suitable for modifying the ^Li abundance, are 
shown as red and blue shaded regions in each panel of Fig. 3. We see that they straddle 
the erstwhile WMAP strips If we had been able to allow a Deuterium abundance as low 
as D/H ~ (1 — 2) X 10“®, the blue shaded region would have been able to resolve the Li 
discrepancy in the context of the CMSSM with gravitino dark matter. The blue region that 
we now regard as excluded by the lower limit on D/H, which is stronger than that used 
in [13], extends to large mi/ 2 . The red shaded region, which is consistent even with this 
limit on D/H, but yields very large ^He/D, is close to the Higgs lower limit on mi /2 for small 
m 3 / 2 , moving to larger mi /2 for larger 1713 / 2 , so as to keep tnsp within the desired range. 

We displayed in Fig. 1 the impact of the improved lower limit on D/H and the new ^He 
constraint on the abundance of an unstable particle, as a function of its lifetime. Interpreting 
this as a constraint on NSP decay into a gravitino, the panels in Fig. 3 show as solid red 
lines the additional restrictions these constraints impose on the (mi/ 2 , mo) planes for different 
values of 7713 / 2 . The effects for small 7713/2 = 10 GeV [in panel (a)] 7713/2 = 100 GeV when 
mo is large [in panel (b)] and 0.27T7 o [in panel (c)] are relatively modest. This is because the 
limit occurs in a region where the NSP is a neutralino, and the relic density varies relatively 
rapidly. Hence a small change in the 7711/2 or mo results in a large change in (x, and the 
old and new bounds are relatively close. However, they do bite in the neighbourhood of 
the shaded ^Li blobs, and have the effect of excluding them entirely. The effects for large 
m 3/2 = 100 GeV and small mo [in panel (b)] and 7773/2 = mo [in panel (d)] are relatively 
large, mainly due to the slow variation of (x near the limit which is characteristic of the f 
NSP region, as can be seen from the relatively wide WMAP strips in this region. Reflecting 
this wide separation between the (old) dashed grey-green lines and the (new) solid red lines, 
we see large red and blue swaths in panel (d), where the ^Li abundance could have been 
reduced, but the ^He/D and/or D/H ratios are unacceptable 

Fig. 4 shows the corresponding ( 7771 / 2 , 7770 ) planes for the choice tan/? = 57, which is 
among the larger values allowed in the context of the constrained MSSM. Although the 
shapes of the allowed regions are rather different from the previous tan/? = 10 case, the 
qualitative conclusions are similar. The ^He constraint again has relatively modest impact. 
However, in all cases the red shaded ^Li regions are excluded by the ^He constraint, and the 
blue regions by the D abundance. We note that, in models with 7773/2 = mo as in Fig. 4(d), 
very little of the ( 7771 / 2 , 7770 ) plane admits gravitino dark matter and, in the viable corner, 
there were no possibilities for depleting ^Li. 

Finally, we consider very constrained models motivated by minimal supergravity (mSUGRA), 
in which not only is 7773/2 = mo, but also the trilinear soft supersymmetry-breaking pa¬ 
rameter A determines the bilinear soft supersymmetry-breaking Higgs-mixing parameter: 

B = A — 7773/2. This is compatible with the values of /i and mx specihed by the electroweak 
vacuum conditions for only one value of tan/? for any given pair of values of ( 7771 / 2 , mo) [17]. 
Such mSUGRA models are then specified by a choice of A = A 7773 / 2 : the panels of Fig. 5 

^This raises the possibility that a discovery of supersymmetry might have an ambiguous interpretation - 
neutralino LSP or gravitino LSP - in the absence of supplementary information. 

^We note in passing that in several panels there are parts of the blue and/or red regions at large mo that 
are also excluded by the WMAP relic density limit. 
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Figure 4: As in Fig. 3, but now for tan/? = 57. 
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Figure 5: ^45 in Fig. 3, but now for very constrained models motivated by mSUGRA. The 
value o/tan/? is fixed by the vacuum conditions, and varies across the {mi/ 2 ,mo) planes with 
values indicated by the steep black contours. These models are specified by the choices of 
A = Am 3 / 2 .' (a) A = (3 — VS), the Polonyi model, and (b) A = 2. 


assume (a) A = (3 — \/3) as found in the simple Polonyi model of supersymmetry breaking 
in mSUGRA [36], and (b) A = 2. For A = (3 — v^), the contours of constant tan/? are 
approximately vertical, and range from about 10 at low mi /2 to about 30 at high values of 
mi/ 2 . In the interesting region of panel (a) where ^Li can be depleted, tan/3 ~ 20 — 30. As 
before, we consider here only regions of the (mi/ 2 ,^ 0 ) planes between the two solid black 
lines: above the higher one, the gravitino is no longer the LSP, and below the lower one 
the lifetime falls below 10^ s. In addition to the constraints discussed earlier, panel (b) also 
displays a small green shaded region at low mi /2 that is excluded by 6 —>• sy decay. 

We see that there is a large difference between the effects of implementing the old and new 
light-element constraints in panel (a) of Fig. 5 for A = (3 — v^), the Polonyi value, whereas 
the effect in panel (b) for A = 2 is smaller. In the Polonyi case, there are large ^Li-friendly 
regions that are excluded by the ^He and D constraints. This reflects the fact that tan (3 is 
relatively small in this case, so the model is qualitatively similar to the 717 . 3/2 = tuq case for 
tan/3 = 10 shown in panel (d) of Fig. 3. On the other hand, tan/3 is typically larger for 
A = 2, and when combined with a smaller stau mass, we find no visible ^Li-friendly region, 
as a result of small yet significant shifts in the values of both C^x and tx- 

6 Conclusions 

In the absence of a convincing astrophysical explanation for the apparent discrepancy be¬ 
tween the observed abundance of ^Li and that calculated on the basis of the baryon-to-photon 
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ratio inferred from CMB observations, it has been natural to explore the possible effects of 
late-decaying massive particles. Indeed, they could suppress the primordial abundance of 
^Li, but at some price. Either the abundance of Deuterium should be very low compared 
with the latest available measurements, and/or the primordial ^He/D ratio must have been 
very high. 

The latest observations of remote cosmological clouds along the lines of sight of high- 
redshift quasars suggest that D/H > 2.2 x 10“^, ruling out much of the parameter space for 
unstable particles that would otherwise have been suitable for diminishing the ^Li abundance 
to agree with observations. A signihcant part of this parameter space would have been 
allowed by the more relaxed limit D/H > 1.3 x 10“^ considered previously [13]. 

The remaining part of the parameter space for unstable particles that is consistent with 
the current lower limit would yield a ^He/D ratio at least an order of magnitude higher 
than the proto-solar value. Since D has been destroyed by stars, reducing its abundance 
from approximately 2.5 to 1.5 xl0“®, while the ^He abundance is thought to have remained 
roughly constant, a primordial ratio of ^He/D > 1 or 2 is unacceptable. This closes the 
remaining loophole for suppressing ^Li without running into conflict with the other light- 
element abundances. 

We have also analyzed the potential bounds imposed by the ®Li abundance. In principle, 
the high growth of its abundance between BBN and the formation of halo stars could also 
have been explained by the decays of the NSP in supersymmetric models, were it not for 
either the low resulting abundance of D/H or the high ratio of ^He/D. 

So what is the interpretation of the apparent discrepancy between the calculations of 
the primordial ^Li abundance and the Spite plateau? We have argued that the origin of 
the discrepancy cannot be the possible existence of unstable particles able to destroy the 
primordial nucleus. This leaves the problem open. Systematic uncertainties in nuclear 
effects such as higher ^Be -|- D reaction rates have been considered [6, 29], but seem unable 
to modify substantially the abundance of ^Li. Stellar mechanisms of depletion may be the 
last resort [26]. Perhaps other new and exciting astrophysical or physical effects will have to 
be considered. 

A signihcant output of this analysis has been the demonstration of the importance of 
the ^He/D constraint on late-decaying massive particles as argued in [14]. Can one make 
the ^He constraint more precise? This would require considering in more detail the cosmic 
evolution of D and ^He. One should allow for the possibility of exotic effects such as large- 
scale destruction of this isotope in primitive structures such as massive Population-HI stars, 
followed by moderate production by normal galactic evolution. However, any such scenario 
should consider simultaneously the cosmic evolution of D. As it is more fragile than ^He, D 
would also be destroyed in any Population-HI stars. As we have seen, the ^He constraint 
sharpens the embarrassment of supersymmetric models with heavy particles whose lifetimes 
are ~ 10^ to 10® s. This interest motivates more detailed studies of the cosmic evolution of 
®He and D, as well as ®Li and ^Li. 

Finally, it would be very useful if the observed abundances of D/H in quasar absorption 
systems were improved. There is currently considerable dispersion in the observed abun¬ 
dances. We note that two such systems in the directions of Q2206-199 at z = 2.0762 with 


14 



D/H = 1.65 ±0.35x10-^ [37] and PKS1937-1009 at z = 3.256, with D/H = 1.6l[]J^x [38] 

have quite low abundances of D, similar to that observed at present in the ISM of our 
Galaxy [39]. If these measurements were to represent the correct D/H abundance in those 
clouds, ^Li depletion and ®Li production by sparticle decay would be a viable option, though 
by solving one problem we would open two new problems. Why is the D/H abundance in 
most of the other absorption systems signihcantly higher, and how can we account for the 
D/H abundances in the solar system, which are also in the range 1.5 - 2.5 xl0“®? 


Acknowledgments 

We thank R. Cyburt, B. Fields, Y. Santoso and V. Spanos for collaborations on related 
topics. The work of K.A.O. and E.V. was supported by the Project ’’CNRS/USA”, and the 
work of K.A.O. was also supported partly by DOE grant DE-FG02-94ER-40823. 


References 

[1] T. P. Walker, G. Steigman, D. N. Schramm, K. A. Olive and K. Kang, Astrophys. J. 
376 (1991) 51; K. A. Olive, G. Steigman, and T. P. Walker, Phys. Rep. 333 (2000) 389; 
B. D. Fields and S. Sarkar, Phys. Rev. D66 (2002) 010001. 

[2] G. Steigman, D.N. Schramm, and J. Gunn, Phys. Lett. B66 (1977) 202; for a recent 
analysis see: R. H. Gyburt, B. D. Fields, K. A. Olive and E. Skillman, arXiv:astro- 
ph/0408033. 

[3] G. L. Bennett et ai, Astrophys. J. Suppl. 148 (2003) 1 [arXiv:astro-ph/0302207]; 
D. N. Spergel et ai, Astrophys. J. Suppl. 148 (2003) 175 [arXiv:astro-ph/0302209]. 

[4] R. H. Gyburt, B. D. Fields and K. A. Olive, Phys. Lett. B 567 (2003) 227 [arXiv:astro- 
ph/0302431]. 

[5] R. H. Gyburt, B. D. Fields and K. A. Olive, New Astron. 6 (1996) 215. 

[6] A. Goc, E. Vangioni-Flam, P. Descouvemont, A. Adahchour and G. Angulo, Astrophys. 
J. 600 (2004) 544 [arXiv:astro-ph/0309480] 

[7] P. Descouvemont, A. Adahchour, G. Angulo, A. Goc and E. Vangioni-Flam, ADNDT 
88 (2004) 203 [arXiv:astro-ph/0407101]. 

[8] A. Guoco, F. locco, G. Mangano, G. Miele, O. Pisanti and P. D. Serpico, Int. J. Mod. 
Phys. A 19 (2004) 4431 [arXiv:astro-ph/0307213]. 

[9] R. H. Gyburt, Phys. Rev. D 70 (2004) 023505 [arXiv:astro-ph/0401091]. 

[10] K. Jedamzik, Phys. Rev. D 70 (2004) 063524 [arXiv:astro-ph/0402344]; K. Jedamzik, 
Phys. Rev. D 70 (2004) 083510 [arXiv:astro-ph/0405583]. 


15 



[11] J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. D 68 (2003) 063504 [arXiv:hep- 
ph/0306024]. 

[12] D. Findley, Astrophys. J. 294 (1985) 1; J. R. Ellis, D. V. Nanopoulos and S. Sarkar, 
Nucl. Phys. B 259 (1985) 175; D. Findley, Phys. Fett. B 171 (1986) 235; M. H. Reno 
and D. Seckel, Phys. Rev. D 37 (1988) 3441; S. Dimopoulos, R. Esmailzadeh, F. J. Hall 
and G. D. Starkman, Nucl. Phys. B 311 (1989) 699; J. Ellis et al, Nucl. Phys. 
B 337 (1992) 399; M. Kawasaki and T. Moroi, Prog. Theor. Phys. 93 (1995) 879 
[arXiv:hep-ph/9403364] and Astrophys. J. 452 (1995) 506 [arXiv:astro-ph/9412055]; 
E. Holtmann, M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 60 (1999) 023506 
[arXiv:hep-ph/9805405]; M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 63 (2001) 
103502 [arXiv:hep-ph/0012279]; K. Kohri, Phys. Rev. D 64 (2001) 043515 [arXiv:astro- 
ph/0103411]. 

[13] R. H. Cyburt, J. R. Ellis, B. D. Fields and K. A. Olive, Phys. Rev. D 67 (2003) 103521 
[arXiv:astro-ph/0211258]. 

[14] M. Kawasaki, K. Kohri and T. Moroi, arXiv:astro-ph/0402490 and arXiv:astro- 
ph/0408426. 

[15] J. F. Feng, S. F Su and F. Takayama, Phys. Rev. D 70 (2004) 063514 [arXiv:hep- 
ph/0404198]; J. F. Feng, S. Su and F. Takayama, Phys. Rev. D 70 (2004) 075019 
[arXiv: hep-ph/0404231]. 

[16] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Fett. B 588 (2004) 7 
[arXiv:hep-ph/0312262]. 

[17] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Fett. B 573 (2003) 162 
[arXiv:hep-ph/0305212]; J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. 
Rev. D 70 (2004) 055005 [arXiv:hep-ph/0405110]. 

[18] K. A. Olive and E. D. Skillman, Astrophys. J. 617 (2004) 29 [arXiv:astro-ph/0405588]. 

[19] D. Kirkman, D. Tytler, N. Suzuki, J. M. O’Meara and D. Fubin, Astrophys. J. Supp. 
149 (2003) 1 [arXiv:astro-ph/0302006]. 

[20] S. G. Ryan, T. G. Beers, K. A. Olive, B. D. Fields, and J. E. Norris Astrophys. J. Fett. 
530 (2000) F57 [ arXiv:astro-ph/9905211]. 

[21] M. Asplund etal. , (2004) in preparation. 

[22] D. F. Fambert, AIP GonF Proc. 743 (2005) 206 [arXiv:astro-ph/0410418]. 

[23] T. M. Bania, R. T. Rood, and D.S. Balser, Nature 415 (2002) 54. 

[24] J. Geiss and G. Gloeckler, Sp. Sci. Rev. 84 (1998) 239; G. Gloeckler and J. Geiss, Sp. 
Sci. Rev. 84 (1998) 275. 


16 



[25] C. P. Deliyannis, P. Demarque and S. D. Kawaler, Astrophys. J. Supp. 73 (1990) 21. 

[26] S. Vauclair and C. Charbonnel, Astrophys. J. 502 (1998) 372; M. H. Pinson- 
neault, T. P. Walker, G. Steigman and V. K. Narayanan, Astrophys. J. 527 (1998) 
180 [arXiv:astro-ph/9803073]; M. H. Pinsonneault, G. Steigman, T. P. Walker and 
V. K. Narayanan, Astrophys. J. 574 (2002) 398 arXiv:astro-ph/0105439; O. Richard, 
G. Michaud and J. Richer, A.A. 431 (2005) 1 [arXiv:astro-ph/0409672]. 

[27] S. G. Ryan, J. E. Norris and T. G. Beers, Astrophys. J. 523 (1999) 654 [ arXiv:astro- 
ph/9903059]. 

[28] J. Melendez and I. Ramirez, Astrophys. J. 615 (2004) L33 [arXiv:astro-ph/0409383]. 

[29] R. H. Gyburt, B. D. Fields and K. A. Olive, Phys. Rev. D 69 (2004) 123519 [arXiv:astro- 
ph/0312629]. 

[30] B. D. Fields and K. A. Olive, New Astronomy 4 (1999) 255; E. Vangioni-Flam, M. 
Gasse, R. Gayrel, J. Audouze, M. Spite and F. Spite, New Astronomy 4 (1999) 245. 

[31] E. Rollinde, E. Vangioni-Flam and K. A. Olive, arXiv:astro-ph/0412426. 

[32] E. Vangioni-Flam, K. A. Olive, B. D. Fields and M. Gasse, Astrophys. J. 585 (2003) 
611 [arXiv:astro-ph/0207583]. 

[33] J. R. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Srednicki, Phys. Lett. B 510 (2001) 
236 [arXiv:hep-ph/0102098]; V. D. Barger and G. Kao, Phys. Lett. B 518 (2001) 117 
[arXiv:hep-ph/0106189]; L. Roszkowski, R. Ruiz de Austri and T. Nihei, JHEP 0108 
(2001) 024 [arXiv:hep-ph/0106334]; A. B. Lahanas and V. G. Spanos, Fur. Phys. J. G 
23 (2002) 185 [arXiv:hep-ph/0106345]; A. Djouadi, M. Drees and J. L. Kneur, JHEP 
0108 (2002) 055 [arXiv:hep-ph/0107316]; U. Ghattopadhyay, A. Gorsetti and P. Nath, 
Phys. Rev. D 66 (2002) 035003 [arXiv:hep-ph/0201001]; H. Baer, G. Balazs, A. Belyaev, 
J. K. Mizukoshi, X. Tata and Y. Wang, JHEP 0207 (2002) 050 [arXiv;hep-ph/0205325]; 
R. Arnowitt and B. Dutta, arXiv:hep-ph/0211417; J. R. Ellis, K. A. Olive, Y. Santoso 
and V. G. Spanos, Phys. Lett. B 573 (2003) 163 [arXiv:hep-ph/0308075]. 

[34] J. R. Ellis, K. A. Olive, Y. Santoso and V. G. Spanos, Phys. Lett. B 565 (2003) 176 
[arXiv:hep-ph/0303043]. 

[35] H. Baer and G. Balazs, JGAP 0305 (2003) 006 [arXiv;hep-ph/0303114]; A. B. La¬ 
hanas and D. V. Nanopoulos, Phys. Lett. B 568 (2003) 55 [arXiv:hep-ph/0303130]; 
U. Ghattopadhyay, A. Gorsetti and P. Nath, Phys. Rev. D 68 (2003) 035005 [arXiv:hep- 
ph/0303201]; G. Munoz, arXiv:hep-ph/0309346; R. Arnowitt, B. Dutta and B. Hu, 
arXiv: hep-ph/0310103. 

[36] J. Polonyi, Budapest preprint KFKI-1977-93 (1977); R. Barbieri, S. Ferrara and G.A. 
Savoy, Phys. Lett. 119B (1982) 343. 


17 



[37] M. Pettini and D. V. Bowen, Astrophys. J. 560 (2001) 41 [arXiv:astro-ph/0104474]. 

[38] N. H. M. Crighton, J. K. Webb, A. Ortiz-Gill and A. Fernandez-Soto, MNRAS 355 
(2004) 1042 [arXiv:astro-ph/0403512]. 

[39] H. W. Moos et ai, Astrophys. J. Suppl. 140 (2002) 3 [arXiv:astro-ph/0112519]. 


18 



